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' Abstract. The dependence of the gluon and ghost propagator in pure SU(3) gauge theory on the 

(N ■ choice of Gribov copies in Landau gauge is studied. Simulations were performed on several lattice 

■ sizes at /3 = 5.8, 6.0 and 6.2. In the infrared region the ghost propagator turns out to depend on the 

D ' choice, while the impact on the gluon propagator is not resolvable. Also the eigenvalue distribution 

, of the Faddeev-Popov operator is sensitive to Gribov copies. 
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Studying non-pertubative features of QCD such as confinement, there are two com- 
mon approaches: lattice gauge theory and Dyson-Schwinger equations. From the latter 
^ I approach there are promising results in recent years [1] about the infrared behavior of 

the gluon D and the ghost propagator G. Denoting by Z the dressing functions of the 
^ ! corresponding propagator, in Landau gauge they can be written as 

|: and G(,^)=5^^. (1) 

(D : 

^ ■ According to [1] in the low-momentum region the dressing functions are proposed to 

behave as Zgi <x (q^)^'^ and Zgh°^{q^)^'^ with a common value K" G (0.5, 1). The in- 
; frared suppression of the gluon propagator and the enhancement of the ghost propagator 

j_| ■ at low-momentum is in agreement with the Zwanziger-Gribov horizon condition [2-4] 

as well as with the Kugo-Ojima confinement criterion [5]. 

Zwanziger [2] has suggested that in the continuum the behavior of both propagators 
in Landau gauge results from restricting the gauge fields to the Gribov region Q., where 
the Faddeev-Popov operator is non-negative. Generically, one gauge orbit has more than 
one intersection (Gribov copies) within the Gribov region Q., but expectation values 
taken over this region are proposed to be equal to those over the fundamental modular 
region A. On a finite lattice, however, this is not expected [2]. In this contribution we 
assess the importance of the Gribov ambiguity on a finite lattice for the SU (3) ghost and 
gluon propagators as well as for the lowest eigenvalues of the Faddeev-Popov operator. 

To study these propagators in Landau gauge using lattice simulation, all thermal- 
ized gauge field configurations {Ux,fi} have to be fixed to this gauge. On the lattice 



Talk presented by A. Sternbeck at Quark Confinement and the Hadron Spectrum VI, Villasimius, 
Sardinia, Italy, September 21-25, 2004. 



90 



90 



2 

[GeV] 



[GeVl 



FIGURE 1. The upper parts show the dressing functions of the ghost Zgf, and gluon propagator Zgi 
measured on best gauge copies as functions of the momentum q (scaled to physical units at )3 = 5.8,6.0 
and 6.2) using various lattice sizes. The lower parts show the ratio (z'*^'^^}/(z'*"^^) determined from the 
first (fc) and best (be) gauge copies. 

the Landau gauge condition is implemented by searching for a gauge transformation 
= SxUx,iig\j^j^ , while keeping 1/^,^ fixed, which maximizes the functional 
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This functional has many different local maxima whose number increases as the lattice 
size increases or the inverse coupling /3 decreases. The different gauge copies corre- 
sponding to those maxima are called Gribov copies, due to its relation to the Gribov 
ambiguity in the continuum [4] . All Gribov copies {^U } belong to the gauge orbit cre- 
ated by U and satisfy the lattice Landau gauge condition dp^^Ax^^ 



with 
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traceless 
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In the literature it is widely accepted that the gluon propagator does not depend on 
the choice of Gribov copy, while an impact on the SU (2) ghost propagator has been 
observed [6-8]. However, in a more recent investigation [9] an influence of Gribov 
copies on the SU (3) gluon propagator has been demonstrated, too. 

Here we report on a combined study of the SU (3) gluon and ghost propagator in 
Landau gauge on the same gauge field configurations generated at /3 = 5.8, 6.0 and 6.2. 
For each configuration we have taken Ncp = 30, 40 and 10 random gauge copies for the 
lattice sizes 16^, 24^ and 32"^, respectively. A subsequent gauge-fixing was carried out 
using standard over-relaxation until mdiXx{dpi^Ax,p,Y < 10^^^ was reached. 

On each first (fc) and each best (be) gauge copy — that with largest functional value 
among Ncp copies — both the ghost and the gluon propagator have been measured. 
The results are shown in Fig. 1. The upper parts show the dressing functions of both 
propagators measured on the best gauge copies as a function of the momentum q scaled 
to energy units. In order to compare to other studies [9, 10] we have used a^^ = 1.53, 
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FIGURE 2. The frequency H{X ) of the lowest eigenvalues X of the Faddeev-Popov operator is shown. 
Full boxes represent the distribution obtained on the best gauge copies, while empty boxes represent those 
on the first gauge copies. 



1.885 and 2.637 GeV for /3 =5.8, 6.0 and 6.2, respectively. Looking at the lower parts 
of this figure it becomes clear that the ghost propagator is affected by the choice of the 
Gribov copy the more the momentum is decreased. The impact on the gluon propagator 
stays inside the statistical error. For further details we refer to [11]. 

In trying to fit to the proposed power laws of the dressing function at lowest momenta, 
mentioned at the beginning, it turns out the lattice sizes used are to small to confirm such 
a behavior. 

We also calculated the eigenvalue distribution of the Faddeev-Popov operator on the 
first and best gauge-fixed configurations as shown in Fig. 2. Looking at this figure it is 
obvious that the distribution H{X) of the lowest lying eigenvalues A on the best gauge 
copies is slightly shifted towards larger eigenvalues compared to that determined on 
arbitrary first gauge copies. Thus better gauge-fixing seems to increase the gap between 
the lowest eigenvalues and the Gribov horizon. 

All simulations were done on the IBM pSeries 690 at HLRN. We thank R. ALkofer for 
discussions and H. Stiiben for contributing parts of the program code. This work has been 
supported by the DFG under contract FOR 465. A. Sternbeck acknowledges support of the 
DFG-funded graduate school GK 271. 
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